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Abstract 

We derive perturbation analytic expressions for the Casimir free energy and entropy between two 
dissimilar ferromagnetic plates which are applicable at arbitrarily low temperature. The dielectric 
properties of metals are described using either the nondissipative plasma model or the Drude 
model taking into account the dissipation of free charge carriers. Both cases of constant and 
frequency-dependent magnetic permeability are considered. It is shown that for ferromagnetic 
metals described by the plasma model the Casimir entropy goes to zero when the temperature 
vanishes, i.e., the Nernst heat theorem is satisfied. For ferromagnetic metals with perfect crystal 
lattices described by the Drude model the Casimir entropy goes to a nonzero constant depending on 
the parameters of a system with vanishing temperature, i.e., the Nernst heat theorem is violated. 
This constant can be positive which is quite different from the earlier investigated case of two 
nonmagnetic metals. 

PACS numbers: 12.20.Ds, 42.50.Lc, 42.50.Nn 
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I. INTRODUCTION 


In the last few years physical phenomena caused by quantum fluctuations of the elec¬ 
tromagnetic field attracted much attention in both fundamental physics and technological 


Q 


applications. Among other fluctuation phenomena, the van der Waals 111 and Casimir 
forces occupy a highly important place because they manifest itself as a macroscopic inter¬ 
action between closely spaced material bodies. The unified theory of the van der Waals and 


Casimir forces based on quantum electrodynamics was developed by Lifshitz 3] sixty years 
ago, but sufficiently precise measurements have been performed only recently (see Refs. M 

die Casimir forces in 

Bose-Einstein condensation jt)], semiconductors 
etc. 


for a review). These measurements opened prospective applications of t 
nanotechnology jl, 8], 
sitions jl3], graphene microstructures 


10 


16], phase tran- 


Although the Lifshitz theory turned out to be very useful for interpretation of the mea¬ 
surement data, the most precise experiments 


19 


231 ] using the configuration of nonmagnetic 


metallic test bodies demonstrated a disagreement with theoretical predictions if the relax¬ 
ation properties of conduction electrons are taken into account in calculations. The same 


measurement data were found in agreement with theoretical predictions of 


he Lifshitz theory 


19H23I]. Simultaneously, 


if the relaxation prop erties of conduction electrons are disregarded 
it was shown 24-26] that for metals with perfect crystal lattices the Casimir entropy found 


in the framework of the Lifshitz theory violates the third law of thermodynamics (the Nernst 
heat theorem) if the relaxation properties of conduction electrons are taken into account. 
With omitted contribution of relaxation properties, the Nernst heat theorem is satisfied. 
The relaxation properties of free electrons are well described by the Drude model. The the¬ 
oretical approach taking these properties into accoimt in calculations of the Casimir force is 
known as the Drude model approach. The plasma of conduction electrons with no dissipa¬ 
tion is described by the plasma model. Usually it is applicable at high frequencies, which 
are much larger than the relaxation frequency. The theoretical approach disregarding the 
relaxation properties of conduction electrons in calculations of the Casimir force is called 
the plasma model approach. 

A disagreement of the measurement data with the Drude model approach and the viola¬ 


tion of the Nernst heat theorem in this approac 
the literature. Specifically, it was underlined [27, 


i are puzzling and created a discussion in 
that for real metals with some fraction 
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of impurities the Casimir entropy jumps to zero at sufficiently low temperature, i.e., the 
Nernst heat theorem is formally restored. This, however, does not provide a satisfactory 
explanation for a puzzle 29]. It was also shown 30] that at large separations, where the 


Casimir force is classical, the plasma model approach violates the Bohr-van Leeuwen theo¬ 
rem, whereas the Drude model approach is consistent with it. Taking into account that at 
short separations below a micrometer, where measurements are most precise, the relative 
difference in theoretical predictions of the Drude and plasma model approaches is equal to 
only a few percent, a more definitive experimental evidence is highly desirable. 

Such an evidence was provided by recent experiments on measuring the gradient of the 


Casimir force between ferromagnetic metals 


31 


33]. Here, the relative difference between 


the predictions of two approaches is either almost zero (Au-Ni system) or has an opposite 
sign with respect to the Au-Au system. These permit to exclude the role of a ny p ossible 


systematic effect that could plague the theory-experiment comparison. In Refs. 31rl33] the 


plasma model approach was again confirmed. Further value of ferromagnetic (or, synony¬ 
mously, magnetic) metals for the Casimir physics is that they provide the possibility to 
perform an experiment where theoretical predictions of the Drude and plasma model ap¬ 


proaches differ not by a few percent but by a factor of 1000 


34 ] (see also Ref. 35]). The 


first data sets of such an experiment are already reported. They indicate conclusively that 
the Drude model approach is excluded, whereas the plasma model approach is in agreement 
with the data [36 ], 

Taking into account the crucial importance of magnetic properties for the resolution of 
a puzzle formulated above, in this paper we derive the analytic expressions for the Casimir 
free energy and entropy in the configuration of two parallel plates at temperature T made of 
dissimilar magnetic metals. All the results below are obtained in the framework of the Drude 
and plasma model approaches in the form of perturbation expansions. In the plasma model 
approach, the used small parameters are the relative penetration depths of the electromag¬ 
netic oscillations into metals multiplied by the square roots of static magnetic permeabilities. 
In the case of the Drude model approach the ratios of the relaxation frequencies to the first 
Matsubara frequency serve as additional small parameters. The derivations are first per¬ 
formed assuming that the magnetic permeabilities are static and then generalized for the 
case of frequency-dependent permeability. 

The obtained expressions are used to investigate the low-temperature behavior of the 
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Casimir entropy for magnetic metals. It is shown that for the plasma model the Casimir 
entropy goes to zero when the temperature vanishes, i.e., the Nernst heat theorem is satis¬ 
fied. For magnetic metals with perfect crystal lattices described by the Drude model, the 
Casimir entropy goes to a nonzero limit depending on the parameters of a system when the 
temperature goes to zero, i.e., the Nernst heat theorem is violated. We prove that at zero 
temperature the Casimir entropy of magnetic metals described by the Drude model can be 
positive. This is different from the case of nonmagnetic Drude metals where the Casimir 
energy at zero temperature is always negative. 

The paper is organized as follows. In Sec. [TT| we derive perturbation expansions for 
the Casimir free energy and entropy calculated using the plasma model approach in the 
case on constant magnetic permeabilities and investigate the limit of zero temperature. In 
Sec.Hnithe same is done when the Drude model approach is used in calculations. Section HVl 
contains generalization of the obtained results for the case of frequency-dependent magnetic 
permeabilities. In Sec. |V] the reader will find out conclusions and discussion. 


II. PERTURBATION EXPANSIONS OF THE CASIMIR FREE ENERGY AND 
ENTROPY IN THE PLASMA MODEL 

We consider the configuration of two (n = 1,2) parallel thick plates (semispaces) at a 
separation a made of dissimilar magnetic metals characterized by the frequency-dependent 
dielectric permittivities e^ n \uj) and magnetic permeabilities fi^ n \u). The Lifshitz formula 
for the Casimir free energy per unit area of plates written in terms of dimensionless variables 
takes the form { 2 , fj 

k D T ,°° r°° 

Jr ( a,T ) = 8 ydy ^2 ln [ 1 ~ r £ ) ( i Ci,y)r£ ) (iCi,y)e~ v ] • (1) 

Here, ks is the Boltzmann constant, T is the temperature and Q are the dimensionless 
Matsubara frequencies connected with the dimensional ones £; = 2-KlksT/% by the relation 
(1 = ii/u c where u> c = c/(2a). The prime in the first sum on the right-hand side of Eq. (jT]) 
means that the term with / = 0 is divided by two. The second sum is over two independent 
polarizations of the electromagnetic held, transverse magnetic (a = TM) and transverse elec¬ 
tric (a = TE). The reflection coefficients in Eq. (pQ) calculated at the imaginary Matsubara 
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frequencies are given by 


ttm = 




4 n) y - \jy 2 + [4 n) f4 n) - ^C? 

4 n) y + v'» 2 + b”V!” ) -i]c, 2 ’ 
dA - yU + dTh - i]g 
dA + \/» 2 + h ( "Vi <n) - i]Cj 2 


( 2 ) 


where = €^(iQu; c ) and /4 ?t) = /r (n ) (i(iu c ). 

In this section we consider the dielectric permittivity of the plasma model which describes 
the nondissipative gas of free electrons 37]. At the imaginary Matsubara frequencies the 
dielectric permittivities of both plates in the framework of the plasma model are 


e, (n) = 1 + 


M' 

CJp 

¥ 


= 1 + 


UJp 

Ci 


(3) 


where Up L) are the plasma frequencies for the metals of the plates and = UJp 1 ' 1 /uj c . 


Now we use the calculation procedure developed in Ref. 


in the case of nonmagnetic 

(n) 


metals. In this and in the next section we assume constant magnetic permeabilities u, ; = 

a. E ,o 


/ 1 q 1 \ Using the Poisson summation formula adapted for the case of even functions 
can be rewritten in the form 


he 


F{a,T) = ir 2 3 XT / ydy I ^Ccos (ItC) F(C,y), 


167r 2 a 3 

where t = T eS /T = he/ (2ak B T) and 


(4) 


1=0 


F(C,y) = i 1 - r *\K,yy { a\<,y)e y ]. 


(5) 


The term of Eq. ([5]) with / = 0 describes the Casimir energy per unit area E(a) at zero 
temperature. The terms with l > 1 represent the thermal correction to it. Then Eq. (j3J) can 
be written as 

E{a,T) = E(a) + A T E{a,T), (6) 


where 


A rHa,T) = 


he 


167r 2 a 3 


/ y d y d c cos (itc) f(c, y)- 


(7) 


i=i 
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Now we consider separation distances between the plates satisfying a condition 


a > A< n > = 

p , ,(«) 

(jJn 


( 8 ) 


where Ap"' ) are the plasma wavelengths for both plates. In this separation region the char¬ 
acteristic frequency u c is much smaller than Up . We also assume that separation distances 
are so large that the following equality is satisfied. 


A = 




47ra 


< 1. 


(9) 


For examp 


metals 


31 


e. for Ni used in experiments on measuring the Casimir force between magnetic 


33] /io ~ HO and \ p /(2i r) ~ 40 nm. Thus, for two Ni plates the inequality ((9]) is 


satisfied for a > 2 pm. Expanding Eq. (J5]) in powers of a small parameter (|9]) , we obtain 


/"2 2 

F( C, y) = 2 hi (1 - e~ y ) + 2 C +V ^ A 


2 e~ y C 4 + y\ 2 


y(e y -i) (l -e~y) y 2 


( 10 ) 


This equation can be substituted in the right-hand side of Eqs. ((4]) and (JT]) and all integrals 
with respect to ( can be calculated explicitly with the result 
rv 


dc cos </tc) = 4 0, w+y‘’WA+y"'(»)A 2 . 


(i), 


| ( 2 )/ 


( 11 ) 


where the functions A[ 0 \ and are given by 
A?\y) =j In (1 - e~ y ) sin {tty) , 


A W (») = 


y ( e y — 1) [ l 3 t 3 


sin (Ity) y cos (Ity) y 2 sin (Ity) 


fit 2 


It 


A?\y) = 


Ae~ y 


y 2 (1 - e~y) . 


12 sin (Ity) 12 y cos (Ity) 


fit 5 


fifi 


6p 2 sin (Ity) 2p 3 cos (Ity) p 4 sin (Ity) 


fit 3 


fit 2 


It 


( 12 ) 


After substitution of Eqs. (11111 and (TT2]l in Eqs. (J4j) and ((7J), the integrals with respect to y 
are also calculated in the form 


poo py 

/ ydy 
lo Jo 

— r (° A + \ r ( 2 )/ Ua 2 


A%) + A (i) MA + 4 2 ) (»)a 2 
= B;“'W+B i (I) WA + B i (2) (i)A 2 , 


(13) 
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where for the functions b[°\ and B\ 2 ^ we hnd 
Bi°\t) = 2 


1 7r cotli (nit) 

¥b~ 


7T 


2l 3 t 3 2 l 2 t 2 sinh 2 (nit) 


B?\t) = -4 

27T 3 cosh (nit) 


n . . 4 7r 2 


^coth(rt)- —+ —— 


It sinh 2 (nIt) 


B?\t) = 2 


7T 


+ 


2vr 4 


l 5 t 5 sinh 2 (nit) 


3 coth (7r/t) sinh 2 (nit) 


— 2 coth 2 (7r/t) 


+ 


n 3 l 3 t 3 

coth (nit) 1 


12 

-27rZi 


—r In (l 
Z 4 f 4 V 


sinh 2 (nit) ' nit 
x 67t 247t 1 

e_ "“) "73^- 


n 2 l 2 t 2 


l 3 t 3 l 3 t 3 e 2nlt — 1 


6 


nl 5 t 5 


Li 2 (e- Mt ) 


Here, Li n (z) is the polylogarithm function. 

As a result, the thermal correction (J7J) takes the following explicit form: 


Ar-Ua.T) = [b^W + B, (1) (0A + B< 2> (t) A 2 ' 

1=1 


(14) 


(15) 


where the temperature-dependent coefficients B[ 0 \ B^ } and B[ 2 ^ are given in Eq. ([14]) . 
Calculating the negative derivative of Eq. (fl5|) with respect to temperature, one obtains the 
explicit expression for the Casimir entropy. 

For our purposes it is desirable to hnd the asymptotic behavior of the thermal correc¬ 
tion (TT5j) at arbitrarily low temperatures T -C T e g. This corresponds to the condition t > 1. 
Keeping only the largest of the exponentially small contributions in Eq. (fT5j) . taking into 
account that 

lim Li n (A) = z (16) 

M-s-o 

and performing all necessary summations, one obtains 


A T B(a,T) 



he f 

Cr(3) 

n 3 

+ 

2n 

87ra 3 \ 

2 1 3 

90 f 4 

t 2& 

+ A 

Cr(3) 

2n 3 

00 1 
+ 

-e~ 

-2tt t 


t 3 

45t 4 

t 



-A 2 

(r(5) 

t 5 

— 167T" 

] e~ M 

! 

5 


(17) 
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where (r(z) is Riemann zeta function. 

Equation (fI7|) gives the possibility to find the asymptotic behavior of the Casimir entropy 

<9 A T F(a, T) 


S(a,T) = 


dT 


(18) 


when T goes to zero. Omitting the exponentially small contributions, from Eqs. (fT71) and (fTHj) 


we arrive at 


cv rp\ k BT 2 (3Ci?(3) 7r 2 

S( “’ T) = IfaV - 45 T 


+ A 


47T 2 

3C^(3) - —r 


-A 




47T 2 


(19) 


where we have introduced the dimensionless temperature 

T 27t 47 raksT 

t = 27r—— = — = —--. 

T eS t tic 


( 20 ) 


As is seen from Eq. (JT9j) , the Casimir entropy goes to zero when the temperature vanishes 
in accordance with the Nernst heat theorem. One can conclude that the Lifshitz theory 
combined with the plasma model provides thermodynamically consistent description of the 
Casimir interaction between magnetic metals (previously this statement was proved for the 
case of nonmagnetic metal plates 24H26]). 


III. PERTURBATION EXPANSIONS OF THE CASIMIR FREE ENERGY AND 
ENTROPY IN THE DRUDE MODEL 


Here, we consider the Casimir free energy (JT]) with reflection coefficients (J2j), as given 
by the Lifshitz theory 2, 3]. However, instead of the dielectric permittivity of the plasma 
model (l3j) . we use the dielectric permittivity of the Drude model at the imaginary Matsubara 
frequencies 




= 1 + 


07 


(n) 


= 1 + 


07 


(n) 


( 21 ) 


&[& + 7 l “>Cr)] Ci[Ci + 7 ( “ ) (r)]' 

In this equation, qW(T) are the relaxation parameters (relaxation frequencies) of the metals 
of plates. The relaxation parameters depend on the temperature and for perfect crystal 
lattices go to zero faster than the first power of T with vanishing temperature [26|, [39]. 
The dimensionless relaxation parameter is defined as yW(T) = y( n )(T) /oj c . As is seen from 
Eq. (I2TT) . at any (i ^ 0 the plasma model ([3]) can be considered as a limiting case of the 


























Drude model (}2Tj) when yhd goes to zero. Generally this statement is, however, incorrect 
because in the limiting case 7 ^ —> 0 the Drude model along the real frequency axis possesses 
a singularity proportional to S(ca) 40 ], 

Now we substitute the dielectric permittivity f|2Tjl in the reflection coefficients (J2J). For 
convenience in calculations, below we supply all quantities found in Sec. HO using the plasma 
model (j3I) with an index p, and the respective quantities found using the Drude model (j 2 Tj) 
with an index D. Specifically, for the reflection coefficients calculated at zero Matsubara 
frequency using the two models, one obtains 


r T > M 1 (°.'!/) = ’tm > (°'») = 1, r°^\ 0 ,y) = 


^ 1 = r fn) 

^ + 1“ " ’ 


<\0,y) = 


(n) 

v 0 y 


'sjrWV + v 2 


( 22 ) 


(n) 1 

Vo y + 


VpWV + V 2 


For the calculation of the Casimir free energy J 7 d(o-, T) using the Drude model, it is useful 
to present it identically as 


Toia, T ) = T p (a, T ) + F D (a, T ) - Jy(a, T) 


(23) 


and to separate the zero-frequency terms of the last two quantities in the following way: 

k B T 


T D (a,T) = T p (a,T) + 


167ra 2 


x 


dy \In [l - r?«( 0 , y)r%g\o, y)e y 


k B T 

8na 2 


x 


In [l - ^te } ( 0, ?/) r TE ) (0, y) e y 

OO /»oo 

/ ydy^i^i 1 - r a { 1 ) «i^y) r a { 2 ) «hy)e~ y ] 
1 _ 1 J Ci 


(24) 


- Ill [1 - r r a m (i 0 , »K (2 >(i<,, y)e s ]} . 

Note that due to Eq. (1221) the TM contributions at zero Matsubara frequency in both models 
cancel each other. 

We expand the products of the Drude reflection coefficients for two plates up to the first 
powers in small parameters ^ n \T) / Ci- Introducing also the notation (3 n = X^/(Ana) <C 1, 
for the TM polarization we obtain 

^?m 1 } (*Cj, y) r ?M (*Ct> y) = r TM(*C j> y) r TM(Kuy) 

7 (1) CO 


y) 


Ci 


7 (2) (F) 

Ci 


y)i 


(25) 
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where the expansion coefficients are given by 


y) — 


PnQy {ho 0 + Pi k 2 ho n) + 2 (y 2 - (?) 


X 


\/PW + Pl(?(yt l) -!) + ho° 

^"tm id(u y)^TM i((h y) 


piqW-q 0^-1)+^ 


(26) 


+ y 2 


In a similar way, for the TE polarization of the electromagnetic held one finds 

^te \ i (i,y)r?E ) (i(i,y) = r TE (i(i,y)iTE (i(i,y) 


1 {1 \T) 

(i 


R?e (i(i,y) 


7 (2) (T) 

(i 


-^TE (*Ci) h)> 


(27) 


where 




Pn^o l) y 


X 


\//3& 2 + - !) + f'o”’ 

r TE > (iCj, */)»'T < B ) (*Ct,K) 


(28) 


« (i4 n) -1)» 2 - CM"’ -1) 


(n) 

ho 


The logarithms containing the products (j25]) and (j27J) can also be expanded in powers of 
the same small parameters 

In [1 - r^ 1) (iCi,y)r^ 2 \iCi,y)e~ y ] 

= In [1 - rP {1 \i(i,y)rP {2 \i( h y)e- y ] 

7 (1) (T) R i a 1 \iC l ,y)e-y (29) 


As a result, Eq. 


Ci l - r £ (1) (i£,, |/)rS (2) (i(i, y)e~ v 

7 (2) (T) i4 2) «i,y)e-« 

0 1 - r£ (1) (®Ci, h)^a (2 ) «j, h)e _y 

can be written in the form 


Fd{&-,T) — T p (a,T) + J-o(a,T) + J 7 7 (a, T), 


(30) 


where Jqj(a, T) is the contribution at zero Matsubara frequency given by 

k B T 


Tj(a,T) 


167ra 2 


'o 


— In 


y 

P(!) 


dy {in 1 - ^(0, h)ffiE 2) (0, y)e y 


1 ~^te (0,h)4 ( E(0,h)e 


=-2/ 


(31) 
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and Jy, is the contribution of all nonzero Matsubara frequencies 


R r( a ) T) = 


k B T ^ 7 W(T) 


87ra 2 


E 

;=i 


Ci 


x 


Ra\i(i,y)e y 


1 y d yJ2\ P {i )( . r , P (2 )r , * „ 
0 a l-T-£ w (iC/,y>a (*Ci,y)e y 

oo ~ (2) 


^Tyf(T) 
871 -a 2 ^ 0 


( 32 ) 


x 


ydy^2 


Ra\i(i,y)e y 


1 - r^ 1 \iC l ,y)r^ 2 \iC l ,y)e~y 
We consider first the contribution to Eq. (130|) at zero frequency. The first integral in 
Eq. (1311) contains r^E ■ Using Eq. (122|) it can be calculated explicitly. The second integral 
in Eq. (l3Tf containing ^ can be expanded in powers of a small parameter A dehned in 
Eq. (J9J), like this was done in Sec. HU Thus Eq. (|3B can be written as 


Ro(a,T) 


k B T( R (3) 

167ra 2 


1 - 


I; I r ( 1 ) r ( 2 ) 
J-U \ r p r p 


Cr(3) 


- 4A + 12A 2 


Then we consider the contribution T-, on the right-hand side of Eq. 
Eq. (l32]h Expanding it in powers of small parameters f3 n) one obtains 


(33) 

dehned in 




k B T 
8 tt a 2 


E 

z=i 


+ 


1 [2) (T) 

0 


7 ( 1 > (r) 

0 

o 

ydy 


ydy 


PiCty/n 0 Ay 0 


y ( e y — 1 ) e y — 1 


AC?V^o 2) Az/vVo 


+ 


. ( 2 ) 


(34) 


y ( e y — 1 ) e y — 1 

Using Eq. (jSJ) and the dehnition of parameters /3 n , Eq. (jM]) can be presented in the form 

(7 W(T) 


R y(a, T) = 


x 


8na 2 


OO r. 

u 

X) 

dy 


U' 


+ w„« 2 EE 

( 1 ) + V PO ( 2 ) 


Uh 


Ci 


+ 


y 


(35) 


Performing summations and integrations in Eq. (135(1 . one arrives at 

k B T e fi Ci?(3) 


Ry( a i T) = 


x 


8n 2 a 2 


In r + 2 + 


„w l (1| (Q , ,/li7 (2) (r) 

P'O ( 1 ) ' V 


7T 


4Ci?(3) 




r 


ca 


( 2 ) 


(36) 
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Now we substitute Eqs. (133]) and 

T D (a,T) = T p {a,T) + 


in Eq. (1301) and obtain 
^b^Cr(3) 


x 


Li s r), V 


( 1 )^( 2 ) 


P i p 


X 


Cr(3) 

ksTesCni 3) 

87 r 2 a 2 

- In t + 2 + 


167ra 2 


4A + 12A 2 




Mo 


7T 


,,( 1 ) 

UJp 


,( 2 )f 


Mo 


, .( 2 ) 

UJp 


-T 


4«3)T <37) 

It is easily seen that the contribution of Jy, to the Casimir entropy goes to zero with 
vanishing temperature. For perfect crystal lattices at temperatures below liquid helium 


temperature it holds y( n )(T) = (26), )39]. Then from Eq. (I35j) we find that 


dJ r ^ / (a,T) _k B T eS ( R ( 3) 


<9T 


47r 2 a 2 


A) 




x T 


In r + 2 + 


cut 


7T 


( 2 ) 


, .( 2 ) 

UJp 


rj-i2 


IT 


[ r 40*(3) 


4Cfi(3) 


7T 


- eff 


r 


(38) 


and 


lim d:F ^ T) = 0 . 
T— >0 <9T 


(39) 


As a result, for the Casimir entropy calculated using the Drude model from Eq. (137|) one 
obtains. 


S D (a,T) = S p (a,T) 
&bCr(3) 


(40) 


167ra 2 


1 - 


T io I r^M 2 ) 
Ll 3 rp rp 


Cr(3) 


- 4A + 12A 2 


dJ r 1 (a : T) 
dT : 


where S p is given in Eq. (fT9l) and dT^/dT in Eq. (1391) . From Eq. (]40|) it is seen that 


S B (a, 0) = lim S B (a,T) 

T—>0 


(41) 


fenCfi(3) 

167T0 2 


1 - 


Til A 1 ) ( 2 ) 
^1.3 1 r p r p 


Cr(3) 


- 4A + 12A 2 


Thus, the Casimir entropy at zero temperature calculated using the Drude model is not 
equal to zero and depends on the parameters of our system, such as the volume (through the 
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separation distance a) and the properties of the plates (through the magnetic permeabilities 
Po"' ) and the plasma frequencies cup 1 ' 1 ). Taking into account that the Casimir entropy is 
the single separation-dependent contribution to the total entropy of the closed 
arrives to the conclusion that in this case the Nernst heat theorem is violated 


svs 


41 


em, one 


42], 


It is interesting to analyze the result (HT)1 in more detail. In the case of one magnetic 

/n\ 

metal we have r/u = 0 and Eq. (1411) takes a more simple from 


Sd(u, 0 ) — 


^bCr{ 3 ) 

167ra 2 


[l - 4A + 12A 2 ] < 0, 


where from Eq. 


A = 






47ra 


(42) 


(43) 


One can see that in this case the Casimir entropy at zero temperature is always negative. 


For two nonmagne 
metals 


fill 


24 


26J. 


ic plates /Iq 11 = 1 and we reobtain the known result for two nonmagnetic 


The most interesting is the case of two magnetic metals. Here, the dependence of So{a, 0) 
on the magnetic permeabilities of the plates leads to unexpected results. For the sake of 
simplicity, we consider similar plates made of magnetic metal with the magnetic permeability 
ho = ho 1 * = ho^ • For two similar plates we have 


ho 


H ho + 1 

2 i 

M (ho + l) 2 


ho + 1 

« 1 -A. 

ho 


(44) 


Taking into account that according to Eq. (1441) 


r 2 n ^ 

p ~ 


4n 

1 —, 

ho 


we find 


Substituting Eq. 


°° r 2n 


Li 3 


71=1 


27T 2 

3ho 


in Eq. (l44jh it is easily seen that under the condition 
.2 

< A 


7r 


(45) 


(46) 


(47) 


6hoCfl(3) 

the entropy at zero temperature is positive. Taking into account the definition of A in 
Eq. (J9]), one arrives from Eq. (|47ll to an equivalent condition 

3A P ho /2 G?(3) 


a < 


7T J 


(48) 
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As an example, for Ni the right-hand side of Eq. (l48|i is equal to approximately 34 /x m. 
Thus, for Ni plates the Casimir entropy at T = 0 is positive over wide range of separations 
from approximately 2 pm (the application condition of our perturbation approach) to 34 fi m. 
This is quite different from the previously investigated case of nonmagnetic metals. 


IV. THE ROLE OF FREQUENCY DEPENDENCE OF MAGNETIC PERME¬ 
ABILITY 


In this section we consider frequency-dependent magnetic permeabilities // n )(i£). It is 


known 


43] that permeability of ferromagnetic metals calculated along the imaginary fre¬ 


quency axis decreases with the increase of £ and at some value £ c specihc for each metal 
abruptly drops to unity. As was noted in Ref. 44j, at room temperature the inequality 
£i 3> £ c holds. Because of this, in all applications of the Lifshitz theory at T = 300 K 
one can put p, at all l > 1 and take the ferromagnetic properties into account only in 


31 


33] in the comparison of the 


the zero-frequency term l = 0. This approach was used 
theoretical predictions with the measurement data. 

The rate of decrease of p^ n '(i£) with increasing £ depends on the value of electric resis¬ 
tance. The lower is the resistance of a ferromagnetic material, the lower is the frequency £ c 
at which /i^(i£) drops to unity 43j. For typical ferromagnetic metals £ c is of order 10 5 Hz. 
Thus, not only at room temperature, but even at relatively low temperature T > 0.001K one 
can put p^(i£) = 1 for all / > 1. This means that for real ferromagnetic metals the thermal 
correction to the Casimir energy calculated using the plasma model can be presented in the 
form 

U °° 

At-F( a, T) = —^ £ [dV) + dV)A. + B, W («)A?' 

1=1 


k B T 
167 T a 2 


V dy In 


1 - 


(i) 

K> y 


k'Wt + s 2 


4A + \! khd 1 ’) 2 + » 2 


X 


( 2 ) 

y 


k 2 ) ( 4 2) ) 2 + » 2 


,-y 


_ l _ i /,,( 2 ) Z ~( 2 )\2 I „,2 

Vo y + V Vo [up ) +y 


( 49 ) 


y dy In 


1 - 


y 


lu . 


(!)>, 


+ y 2 y 


Ul' 


( 2 )> 


+ y 2 


,-v 


y+ yj+ y+ \J(u ^) 2 + y 2 
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where 


Ai 


Ap 1} + A p 

4n a 


( 2 ) 


<C 1. 


(50) 


The sum on the right-hand side of Eq. (|49|) is the thermal correction (fT5]l , where the mag¬ 
netic properties are omitted [this is seen from the replacement of A defined in Eq. (J9]) with Ai 
defined in Eq. (1501) ]. The following two integrals add the contribution of the zero-frequency 
term with included magnetic properties to the thermal correction (the first one) and sub¬ 
tract the contribution of the same term with omitted magnetic properties. [Note that the 
replacement of // n )(i£) with unity in the frequency region [0,£ c ] leads to only a negligibly 
small influence on the value of E(a).] In Eq. (T49|) we have also taken into account that 
^tm (0) y) = 1 for both magnetic and nonmagnetic materials. As a result, only a difference 
in the values of r^(0,y) for // = and /i = 1 contributes to A T J r {a 1 T). 

We expand the first and second integrals in the right-hand side of Eq. in powers of 
small parameters A and Ai, respectively, and perform integrations with respect to y. The 
result is 


, oo 

T) = TfSv 51 K’w + B ?MAi + 4 

1=1 (51) 

+ w (A_Al)|1 ^ 2(A+Al)1 ' 

For the case of similar plates we have 

A ± Ai = 1) ——. (52) 

Ana 

Note that alternatively we can do not expand the integrals on the right-hand side of Eq. (TI9|) 
in powers of small parameters, but calculate them numerically. In this case the application 
region of Eq. (T49|) is determined by the condition Ai <C 1. For example for Ni this leads 
to a > 200 nm [i.e. Eq. (l49j) is applicable starting from an order of magnitude smaller 
separations than Eq. f[5T]) ]. 

One can use Eqs. (HHj) and (I5T|) to calculate the thermal correction to the Casimir energy at 
T > 0.001K. From Eq. (1491) and Eq. (fl4l) it is easy also to obtain the respective expression for 
the thermal correction to Casimir pressure between two parallel plates made of ferromagnetic 
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metals 


A T P(a,T) = —^-A T F(a,t) 
oa 


he 


Ai 


7 r 3 coth (7r/t) 
87r 2 a 4 [ ( lt) A It sink 2 (nIt) 

1 


E 


7 r 


It sinh 2 (7r/t) _ (nit) 


sinh (nit) cosh (nit) 


+ 4 cotli (nIt) + 2nlt — 6nlt coth 2 (nIt) -f-— 

nit 


3A 2 


n 


[—Anlt + 5 (nIt) 2 coth (nIt) ] 


It sinh 2 (nIt) 

12nlt coth 2 (nit) — 8 (nIt) 2 coth 3 (nit) — 4 coth (nit) 


nksT 

16a 3 


(A - AO 


1--(A + A0 


(53) 


We remind that it is not possible to consider the limiting case T —> 0 in Eqs. ([49]) . (]5T]) 
and (153|) because these equations are obtained under a condition that the temperature is 
larger than some hxed (small) value. In order to investigate the role of the frequency 
dependence of //(*£) at T —^ 0 it is convenient to use the Abcl-Plana formula similar to 
Ref. 45] (see also Ref. Q) 


E'*w 


1=0 


/** , , , $(it) -$(-it) 

/ Ht)dt + i dt -- 

'o Jo 


e 2nt - 1 


Now we choose 


where F((, y) is c 


4(0= / ydyF( C,</), 


(54) 


(55) 


in the form 


45] 


ehned in Eq. (15j) . and obtain the thermal correction to the Casimir energy 

&(irt) — $(— irt) 


a -r-/ \ ^cr r i 


e 2nt - 1 


(56) 


The dependence of the magnetic permeability [J,(i£) on the frequency is described by the 
Debye formula [43] 

//A - 1 l,n - 1 

(57) 


/•>\ T . A^o — 1 , , ho - 1 

h«) = 1 + -—— = 1 + 


1 + 77- 1 + 3 SmC 

where ae m = 0 J c /co m and uJ m is some characteristic frequency which is different for different 
materials. For simplicity we consider the case of two similar magnetic metals and restrict 
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ourselves by the first order terms in A. Then from Eq. (1550 we obtain 


*(0 = 2 

- 2 


ydy 
(/i 0 - 1) aa 

\fy o 


In (1 - e ~ y ) + 


-CA 


y {e y - 1) 

r °o C 2 + y 2 


( 58 ) 


'C 


e y — 1 


dy. 


The contribution of the first integral on the right-hand side of Eq. (1551) to A T J r (a,T) was 


calculated in Ref. 


45]. Because of this, now we consider only the second one 

°°C 2 + y 2 


MO = -2 (Ao iR- CA 

\/ho 


'c 


e y — 1 


dy. 


(59) 


In the lowest order to the small parameter £ = tr, we have 

r e + y 2 

Jc & - 1 

Then, in the lowest order in r one obtains 


dy = C«(3) + 0(C). 


(60) 


$i (irt) — = —4i 


■ (ho 1) ®m^(ii(3) 


a/ ho 


rt. 


(61) 


Substituting this in Eq. ([56]) together with the result of Ref. 45] for the first integral on the 
right-hand side of Eq. (]58|) . we arrive at 


A rHa,T) = 


he f (n(3) _ 7r 3 

8 vra 3 \ 2 1 3 90 1 4 

Cr{3) (ho - 1) 3s m A 


+ A 


Cfl(3) 2 tt 3 

t 3 45 1 4 


(62) 


6tt ,/ho t 2 

This result is in agreement with Eq. ra. but contains an additional term due to the 
frequency dependence of /i. 

Thus, for magnetic metals, an account of the frequency dependence of /i at small T gives 
rise to the second order in T correction in the free energy. As a result, the entropy (1T91) 
acquires a correction 

( Hr, A 

(63) 


m _ k B Ch(3) (mo - 1) aSmA 

l\o\CL , I J — ;—-—— - T. 


247T 3 a 2 y / /io 

As is seen in Eq. (1631) . this correction goes to zero when T goes to zero, i.e., for the plasma 
model the Nernst heat theorem is preserved even with account of the frequenc y d ependence 
of /i. This is in analogy to the case of magneto dielectrics investigated in Ref. 46|. 

For the Drude model, an account of the frequency dependence of /i does not change 
the fact that for metals with perfect crystal lattices the Nernst heat theorem is violated. 


17 























This follows from Eq. (l30]h where the violation of the Nernst theorem originates from the 
contribution of the zero frequency terms entering J r o(a,T). For Tpia, T) on the right hand 
side of Eq. (j3U|) the Nernst theorem is satisfied for the frequency-dependent /i and J-o(a, T ) 
does not depend on the presence of the frequency dependence. Finally, the term T y (a, T ) 
acquires a correction of higher order in T to the terms written in Eq. . This proves that 
the frequency dependence of /i does not change our conclusions concerning a consistency 
of the plasma and Drude model approaches with thermodynamics for the case of magnetic 
metals. 


V. CONCLUSIONS AND DISCUSSION 


In the foregoing, we have investigated analytic behavior of the Casimir free energy and 
entropy between two parallel plates made of dissimilar ferromagnetic metals. In so doing, the 
dielectric properties of metals were described either by the nondissipative plasma model or 
by the Drude model taking into account the dissipation of free charge carriers at arbitrarily 
low temperature. 

Using the perturbation expansions of the Casimir free energy in small parameters it 
was shown that the Lifshitz theory combined with the plasma model satisfies the Nernst 
heat theorem, i.e., the Casimir entropy goes to zero when the temperature vanishes. Quite 
differently, for ferromagnetic metals with perfect crystal lattices described by the Drude 
model it was shown that in the limit of zero temperature the Casimir entropy goes to a 
nonzero limit depending on the parameters of the system, i.e., the Nernst heat theorem is 
violated. Both constant and frequency dependent magnetic permeabilities are considered. 

As was noted in Sec. [U violation of the Nernst heat theorem for the Casimir ent ropy 


between nonmagnetic metals with perfect crystal lattices is a known effect 


flEi, 


24 


26]. 


The distinctive feature of magnetic metals found in this paper is that the Casimir entropy 
at zero temperature depends on the static magnetic permeabilities of plate metals and can 
be positive (recall that for nonmagnetic metals with perfect crystal lattices described by 
the Drude model the Casimir entropy at T = 0 is always negative). This establishes a link 


between ferromagnetic metals described by the Drude model and dielectrics wit 


account dc conductivity of plate materials. In the latter case it is known flfl 


l taken into 


47 


48 1 that 


the Lifshitz theory violates the Nernst heat theorem and the Casimir entropy at T = 0 is 
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positive. For magnetic dielectrics this positive quantity remains independent on the magnetic 


properties 


46|. This makes unique the case of ferromagnetic metals considered here. Note 


also that we have considered the configuration of two parallel plates, but a violation of the 
Nernst heat theorem for the Casimir entropy calculated using the Drude model holds also 
for other configurations. Thus, for a sphere above a plate made of nonmagnetic metals this 
was demonstrated in Ref. 49], 


It is pertinent to briefly discuss the physical meaning of the Casimir entropy and its 
sign. For this purpose we remind that the Casimir free energy per unit area of plates 
(P) is derived by subtraction of the free energy for infinitely separated (uncoupled) plates 


0 . 


from the nonrenormalized free energy of plates separated by a distance a [21- Then the 
Casimir entropy calculated by Eq. (fT8j) also represents respective difference and has the 
meaning of an entropy of the fluctuating field. It characterizes an interaction between 
the plates, but considers only a minor fraction of the entropy of a closed system which 
includes also much larger entropies of the plates. As a result, at not too low temperature 
the total entropy is always positive irrespective of whether the Casimir entropy is positive 
or negative. An important point, however, is that only the Casimir entropy depends on 
the separation distance, whereas the entropies of the plates do not depend on separation. 
Because of this, the total entropy at T = 0 is separation-dependent, i.e., the Nernst heat 
theorem is violated. In fact the Casimir entropy per unit area of two parallel plates is not 
an immediately measured quantity. It can be experimentally found indirectly by means of 
numerical differentiation from the force between a sphere and a plate measured as a function 
of temperature (in the Derjaguin approximation 2] the latter quantity is proportional to 
the Casimir free energy per unit area of two parallel plates). 


As was mentioned in Sec. HI the most precise experiments using nonmagnetic metals 19 
2.3] and all the experiments using magnetic metals 31 h33 


are in agreement with the 
plasma model approach and exclude the Drude model approach to calculation of the Casimir 

he measurement data of most precise experiments using dielectric 


force. In a similar 

way, t 

re n 

test bodies [9, 

15, 

16, 

50 

-52] 


504521 are in agreement with theoretical predictions of the Lifshitz 


theory only if the dc conductivity of boundary materials is omitted in calculations. In our 
opinion, it cannot be accidental that in so many experiments on measuring the Casimir 
interaction the data are in agreement with thermodynamically consistent theoretical ap¬ 
proach and exclude the approaches where the Nernst heat theorem is violated. Thus, the 


19 




























problem of a proper account of the relaxation properties of free charge carriers in metals and 
the dc conductivity in dielectrics when calculating the Casimir interaction invites further 
investigation. 
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